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We study the parameter space of cold dark matter axions in two cosmological 
scenarios with non-standard thermal histories before Big Bang nucleosynthesis: the 
Low Temperature Reheating (LTR) cosmology and the kination cosmology. If the 
Peccei-Quinn symmetry breaks during inflation, we find more allowed parameter 
space in the LTR cosmology than in the standard cosmology and less in the kination 
cosmology. On the contrary, if the Peccei-Quinn symmetry breaks after inflation, 
the Peccei-Quinn scale is orders of magnitude higher than standard in the LTR 
cosmology and lower in the kination cosmology. We show that the axion velocity 
dispersion may be used to distinguish some of these non-standard cosmologies. Thus, 
axion cold dark matter may be a good probe of the history of the Universe before 
Big Bang nucleosynthesis. 

PACS numbers: 14.80.Va, 95.35.+d, 98.80.Cq 

The standard cosmological model has been tested up to a temperature T ~ 1 MeV, or 
down to times as short as ~ 1 s, when Big Bang nucleosynthesis (BBN) occurred. The success 
of the BBN theory is due to its great precision in predicting the primordial abundance of 
light elements D, ^He and ^Li. For the success of BBN, the Universe must be radiation- 
dominated at temperatures T ^ 4 MeV 0]. However, due to lack of data prior to BBN, the 
history of the Universe in the pre-BBN epoch T ^ 4 MeV is only indirectly inferred. 

In the standard cosmology, radiation has been dominating the energy density of the 
Universe before BBN since the very early time at which inflation ended. How radiation was 
produced at the end of inflation from a state of negligible temperature is still a topic of active 
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research: models include the decay of the inflaton field and parametric resonance 

In alternative cosmological models, inflation may have ended at times close to BBN |5|], 
or there might have been a period after inflation in which the dominant energy density was 
not in radiation but in some other exotic form, like the energy density of a scalar fleld I™, 

on 

or still more there could have been an injection of entropy into the radiation p, |8|. 

A good probe of the history of the Universe before ~ Is is a relic particle that has 
survived from that period. Indeed, in the standard cosmology, dark matter relics like axions 
[qI. [lol. ll| and Weakly Interacting Massive Particles (WIMPs) (see 12|) are produced when 



the Universe was ~ 1/is or between ~ 10ns and ~ lOps old (corresponding respectively 
to the age of the Universe at the QCD phase transition or at the WIMP freeze-out for 
WIMP masses between 10 GeV and 1 TeV). These dark matter relics are therefore excellent 
candidates to test the cosmological history of the Universe at very early epochs is!]. 

In this article we study cold dark matter (CDM) axions as probes of the pre-BBN 

epoch, assuming that axions provide the totality of the CDM observed. 



CDM, 



with ll 



QcDuh^ = 0.1131 ±0.0034. 



(1) 



(2) 



Here Qa and ^cdm are the axion and CDM energy densities in units of the critical density 
Pcrit, and h is the Hubble constant in units of 100 km s~^ Mpc~^. 
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In the LTR cosmology the 



expansion of the Universe after inflation is driven by a massive scalar fleld that decays 
and reheats the Universe. This stage lasts down to a temperature Trh, after which standard 
radiation-dominated cosmology applies. In the kination cosmology, the energy content of the 
Universe is dominated by the kinetic energy of a scalar fleld which evolves without entropy 
release down to a temperature Tkm, after which standard cosmology begins. 

We flnd that with respect to the standard cosmology, if the Peccei-Quinn symmetry in 
the axion theory breaks during inflation, the allowed axion parameter space is enlarged in 
the LTR cosmology and restricted in the kination cosmology. Instead, if the Peccei-Quinn 
symmetry breaks after inflation, the mass of cold dark matter axions in the LTR cosmology is 



3 



orders of magnitude smaller than its standard-cosmology value, and it is orders of magnitude 
greater in the kination cosmology. 

Past work on axions in non-standard cosmologies has examined^ the parameter space of 
hot dark matter axions in the LTR and kination cosmologies [40;], and the cosmological 
bound Qah"^ < 1 for cold (i.e. non-thermal) axions in the LTR cosmology assuming that the 
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Peccei-Quinn symmetry breaks after inflation 
as 100% of the cold dark matter, allows the PQ symmetry to break after or during inflation, 
and includes anharmonicities in the axion potential. 

In Section [U we review the two most important mechanisms of axion production, the 
misalignment mechanism and axionic string decays, focusing on how to extend the standard 
results for these mechanisms to any non-standard cosmology. In Section [Tll we review the 
present axion energy density Q^J^'^h'^ in the standard cosmology. In Section [XTTl we derive the 
density Q^^^h"^ and present our results for the axion parameter space for the LTR cosmology. 
In Section IIVI we do the same for the kination cosmology. In Section |V] we compare our 
results to previous work, discuss how our results may be modified by different choices of 
axionic string parameters, and suggest a possibility of distinguishing non-standard pre-BBN 
cosmologies observationally. Our conclusions are presented in Section IVII 

I. PRODUCTION OF AXION DARK MATTER 



The axion is a hypothetical pseudoscalar particle first introduced in the QCD sector of 



ll|. Axion properties depend 
43|). They are: the scale /„ 



the Standard Model to solve the strong CP problem 
on several parameters (see e.g. |4l|] and the recent reviews 42|, 
at which the Peccei-Quinn (PQ) symmetry breaks, the PQ color anomaly A^, the number 
of degenerate QCD vacua A^^^, and the couplings of axions to quarks, leptons, and photons. 
Axion models include the KSVZ or hadronic model 4J], in which there is no tree- level 
coupling between the axion and the standard model quarks and leptons, and the DFSZ 



model 



451, in which the axion couples to standard model quarks and might couple to 



leptons. In both models, the axion-photon coupling is non-zero. Unless otherwise specified, 
our results apply to both KSVZ and DFSZ models (in both of which N = Nh).^ 
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accomplished for specific conditions on the parameters fa and N, on the expansion rate of 
the Universe Hj at the end of inflation, and on the mechanism of axion production. 

In this section, we review the two most important mechanisms to produce axion dark 



matter: vacuum realignment 



and string decays 



471, 



481. We focus on the modifications 



needed to go from the standard cosmology to non-standard cosmologies. 



A. Axions from vacuum realignment 

We review the production of cold axions by the vacuum realignment mechanism, following 



the conventions in Ref. [6j]. The formulas in this section depend on the details of the 
cosmology before BBN only through the dependence of the Hubble expansion rate H{T) on 
the temperature T. Different pre-BBN cosmologies differ in the choice of H{T). 

When the temperature of the Universe falls below T ~ /„, the Peccei-Quinn symmetry 
breaks and the axion field a{x) originates. The equation of motion for the misalignment 
angle 9{x) = a{x)/ fa is 



9 + 3H{T) e+ — v'ie) = 0. (3) 



1 

J a 

Here a dot indicates a derivative with respect to time; V'{6) is the derivative with respect 
to 6 of the axion potential 

V{e) = ml{T)f'A^-^ose). (4) 



The axion mass depends on temperature as 



68] 



1 T<A, 

^a(T) ={ ^ (5) 

bma{^) T>A, 

where we choose = 1, A = 200 MeV, b = 0.018 [57]] and the zero-temperature axion mass 

ma = TUaiT = 0) = 6.2 /ieV /.^V (6) 

We use the notation fa,i2 = GeV. For T ^ A, we can neglect rn\{T) and so V'{6) 

in Eq. ([3]). A solution to Eq. ([3]) in this case is then 9 = 9i = const, where 9i is the initial 
value of the misalignment angle. When T ~ A, coherent oscillations of the axion field begin. 
The temperature Ti at which oscillations being is defined by the condition 

m{T^)=ma{T^). (7) 
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The axion number density at Ti can be written as [49|, |57|, |58|, |59|, |6J] 

n,{T,) = l-ma{T^)f,x{Olfm. 



(8) 



Here angle brackets indicate an average over the possible values of 6f within a horizon 
volume, X = 1-44 is a fudge factor, and f{Oi) is a corrective function accounting for anhar- 
monicities in V{di) \^\. (In our work on cold dark matter axions in standard cosmology 
6J], we used the symbol ( for the parameter x, but here ( is used for another quantity, see 
Section lib.) 

A complete study of the axion parameter space requires taking into account the function 
f{Oi) in the whole rarige — tt < 9i < h 
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various papers 
is [641 



46 
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661]. The function f{Oi) has been studied in 



|69|, |70|. An analytic function interpolating between numerical results 

7/6 



In 



(9) 



The properties of the cosmological axion differ whether the axion field is present during 
inflation (/„ > Hj/2tc) or not (/„ < Hi/2tt). 

For fa < Hi/2'K (Scenario I), the initial misalignment angle 9i acquires different values 
within the same Hubble horizon. Using Eq. ^ we obtain {OffiOi)) = 8.77 [M], where the 
average is taken over the possible values of 6i. 

If fa > Hj/2tt (Scenario II), the PQ symmetry breaks during inflation and it is never 
restored. The initial misalignment angle di takes only one value within a Hubble horizon. 
As noticed by Linde 7l|], to a relatively small value of 6i (~ 10~^) there corresponds a large 
value of fa (~ 10^^ GeV). This region of parameter space can thus accommodate values of fa 
of the order of the GUT scale (~ 10^^ GeV) or above, however at the expense of tuning the 
initial misalignment angle to a small value which was the problem that the PQ mechanism 
tries to solve. 

In Scenario II the mean {6i) = 9i has a unique value within one Hubble volume, so 



{oifm = {o^+^i)f{o. 



Here the variance 



Hi 

2nfa 



(10) 



(11) 



is due to the axion quantum fluctuations present at the end of inflation 
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In Scenario II, axion isocurvature fluctuations are constrained by a combination of Cosmic 
Microwave Background (CMB), Baryon Acoustic Oscillations (BAO) and Supernovae (SN) 



measurements [l5| as 

i// < 4.168 X 10-^ (12) 

This constraint is obtained from an upper limit on the ratio of isocurvature and adiabatic 
power spectra at the wave number = 0.002 Mpc~^. Notice that this bound applies 
independently of the cosmological model between inflation and BBN because modes with 
wave number ko are still well outside the horizon at the epoch of recombination. 

A further bound on Hj comes from the non-detection of primordial gravitational waves. 
The combined CMB+BAO+SN measurements set an upper limit on Hj given by [l^ 

Hj < 6.29 X 10^^ GeV. (13) 

In Scenario I (/„ < Hi/2n), this bound on Hj implies an upper bound on fa equal to 

fa < 1.00 X 10^^ GeV, for fa < Hi/2-k. (14) 

Astrophysical considerations on the cooling time of white dwarfs yield the bound [3] 

/a > 4 X 10^ GeV, (15) 

valid for KS VZ axions. A similar bound from supernovae applies to other axion models [3] • 
The bounds in Eqs. ( fT3l) and ( fTSl) apply independently of the pre-BBN history of the 
Universe. 



B. Axions from string decays 

The spontaneous breaking of the Peccei-Quinn symmetry leads to the formation of topo- 
logical defects such as axionic strings {74]. In Scenario I {fa < Hi/2tt), the PQ symmetry 
breaks after inflation and cold axions produced by the decay of topological strings in the 
early Universe contribute a large fraction of the axion energy density. Instead, in Scenario II 
{fa > Hi /27r), topological defects are washed out by inflation so axions from axionic strings 
are not present. 



The present axion energy density produced by axionic strings p^*''(To) is proportional 



the present axion energy density produced by the misalignment mechanism p™'^(To) 
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to 
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75|. The ratio a between p^*'"(To) and p™'*^(To) can be put in the form 



a 



Here, following the notation in Ref. 



c 



(16) 



52j (however, we use C, for their parameter x to avoid 



confusion with our x in Eq. ([8])): A''^ is the number of degenerate QCD vacua, f is the factor 
by which the axion comoving number density increases due to string decays, averaged over 
all possible processes that convert strings to axions, ^ is a constant factor depending on the 
string network model, and C accounts for the uncertainties in the low-energy cutoff of the 
radiated axion field. 

In the standard cosmology, the numerical values of these parameters have been discussed 



extensively both theoretical 
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55, 



77|, 



ly and via numerical simulations of string networks 47 



48 



5], 



781 ]. However, there is still disagreement about the numerical 
values of f and ^ in the standard cosmology, f^^^ and 

In the following we discuss each parameter in Eq. f|T6l) separately in standard cosmology 
and we extend previous theoretical results to obtain the values of the parameters in a generic 
non-standard cosmology. 

Parameter Nd- In the standard cosmology, it is usually assumed that A'd = 1, because 



for > 1 a domain wall problem may arise 



761]. In modified cosmologies we take the same 



value Nd = 1 as in the standard cosmology, because describes a property of the axion 
field. 

Parameter f. The value of f depends on the details of the axionic string relaxation 
toward lower energy configurations and on the energy spectrum of the radiated axions (53 ]. 
In the standard cosmology it is (see Ref. 521 ]) 



^std 



ln(ti/5), for a slow-oscillating string, 
0.8, for a fast-oscillating string. 



(17) 



Here, ti is the time at which the axion field starts to oscillate and 5 is the string core size 



In Eq. ( 1171) . the first line corresponds to the string emission mode 
while the second line corresponds to the model in Refs. 



48 



in Refs. 



47 



50, 



53 



55, 



The time ti can be expressed in terms of the corresponding temperature Ti using the 
relation 

kt = ^^^^ 
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valid in the standard cosmology. In the illustrative case S = (10^^ GeV) ^ and Ti = 1 GeV, 
one finds f®**^ ~ 70, which is approximately the same value found in Refs. 47|, |50|, l53 



mm 



To extend Eg . ffTTI) to non-standard cosmologies, we repeat its standard-cosmology deriva- 
481 . 1521] but change the relation between time and Hubble parameter in Eq. ( |T8|l 



tion in Refs. 

to that appropriate for a non-standard cosmology. We consider a generic dependence of the 
scale factor a{t) on time t, 

a{t) oc t^, (19) 

where /5 is a constant that depends on the details of the modified cosmology. For example, 
(3 = 2/3 for the matter-dominated Universe and the LTR cosmology, (3 = 1/2 for the 
radiation-dominated Universe, and (3 = 1/3 for the kination cosmology. Eq. ( |T8l) then 
becomes 

m - ^ = T- (20) 

Harari and Sikivie 



the equations 



and 



48| derive the axion number density from string decays n^a'^{t) from 
rfnf (t) 1 dpr 



dt u{t) dt 
dpr dp, 



m{t)nf{t) 



- 2H{t)p,. 



(21) 



(22) 



dt dt 

Here, pr is the energy density of the radiated axions, uj{t) is the average energy of axions 
radiated in string decay processes at time t 
given by 

Ps — ^ J.9 • V^"^J 



521 ]. and ps is the energy density in strings, 
^Nji^flHt/S) 



From Eqs. (EI]), ([22]) and ([231) we obtain 

eiV|/2 27r(l-/3) '"^ 



nfit,) 



,3/3 



dt ln(t/(5) 



(24) 



where tpq <^ti is the time at which the PQ phase transition occurs. The formula to obtain 
the parameter f follows from Eq. (2.13) in Ref. 52l ]. 



nf{t) 



which, evaluated at time ti, gives 



27r(l-/5) 



,3/3-1 



dt \n{t/5) 



(25) 



(26) 
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The function u!{t) depends on the model for the energy spectrum of the emitted axions. 



For slow-oscillating strings, Davis 



471 ] argues that the energy spectrum of the axions radiated 



at time t is peaked around 27r/t, and finds uj{t) = 2ir/t. Using this expression of u{t) in 
Eq. ( l26l) gives 

i4ln(ti/5), for 1/3, 



f ln(ti/tpQ)ln(ti/5), for /3 = 1/3. 



(27) 



In particular, for the standard cosmology (3=1/2 and r*^**^ = ln{ti/6), as in the first line of 
Eq. (E 



For fast-oscillating strings, Harari and Sikivie 48|] argue that the energy spectrum of the 
radiated axions is broad, with a low-energy cutoff at energy ir/ti and a high-energy cutoff 
at energy tt/6. They find uj(t) = (27r/t) ln(t/5). Eq. (!26l) with this expression of uj(t) leads 



to f*^**^ = 1 for the standard cosmology. Numerical simulations [52, 15^, |75[ favor a slightly 
smaller value of f^^*^, namely 0.8 as quoted in Eq. (|T7j) . Therefore, we decided to multiply 
Eq. (!26|) by 0.8. Hence, for the fast-oscillating strings. 



0.8 



1-/3 
3/3-1' 



for p 7^ 1/3, 



0.8| ln(ti/tpQ), for /5 = 1/3. 



(28) 



In the particular case of the LTR cosmology (/? = 2/3), Eq (!28|) gives r^"^^ = 0.27, while in 
the kination cosmology (/? = 1/3), Eq. (EHl) gives f"^'"" = 0.53 ln(ti/tpQ). 

In presenting our results, we use the value of f for fast-oscillating strings, Eq. ( l28l) . We 
discuss the alternative choice of Eq. ( 1271) in Section |Vl 

Parameter ^. The value of in the standard cosmology has been discussed in the 



literature and different authors quote different results 



for an evolving string network in Ref. 
75| give f ~ 1. 



52 



55 



75 



781 ] . Numerical simulations 



78] yield C,^^'^ ~ 13, while simulations in Refs. 



52^, 



55, 



The value of ^ changes in a modified cosmological scenario. The authors in Ref. [3] 
outline a method to estimate ^ in a modified cosmology in which the Universe is matter- 
dominated from the value ^'^^'^ for a radiation-dominated Universe. Here we generalize the 
results of Ref. [3] to a generic cosmology with arbitrary /5, following their method. We 
define the characteristic length L for an axionic string of energy density p and tension n 
per unit length through p = p/L"^. The parameter ^ appears in the time dependence of the 



10 

string length L as L = t/y/^. The method of Ref. (79| consists in computing ^ as 

^ = iloHt)', (29) 
where 70 is the fixed-point value of Eq. (14) in Ref. [3], 

^ = "p + |2^ + 3b-l|, (30) 

Here 7 = {H{t)L)-\ and c> is a constant determined by the value of ^ for a radiation- 
dominated universe. 

For a generic jS, we have 

HI 

and the fixed-point of Eq. fl30l) follows from setting its right-hand side to zero, 

2-3/j+ v/(4c + 9)/32-12/3 + 4 
= W~c • 

Then from Eq. (l29l) and H = f3/t, we find 

e = ^ (2 - 3/? + V(4c + 9)/52-12/3 + 4)'. (33) 

The constant c is fixed from the requirement that ^ = ^^^'^ in a radiation-dominated Universe, 
P = 1/2. This gives 

c-i^. (34) 

This is to be substituted in Eq. ( l33l) to find the value of ^ in the non-standard cosmology. 
For f = 1, we find c = 3/4 and 

e = ^ (2 - 3/5 + 2^3/32-3/5 + 1)' (for e^-i = 1). (35) 

Then, for the LTR cosmology, in which /5 = 2/3, we find = 16/27 = 0.5926; for the 
kination cosmology, in which /5 = 1/3, we find e'^'" = 4(7 + 4^3) /27 = 2.0634. 
The choice ^^^'^ = 13 is discussed in Section IVl 

Parameter C. In the standard cosmology, C ~ 1 [52]. To find C in a generic cosmology 



we use the fact that on dimensional grounds, ( is of order [52], so that to a change 
A,^ there corresponds a change AC/2. However, the theoretical uncertainty on C is around 
50% |5^, higher than the difference A^ due to the change in the cosmology used. We thus 
consider ( constant and equal in all cosmological scenarios. 
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II. AXION CDM IN THE STANDARD COSMOLOGY 

In this section we review the derivation of the present cold axions density, assuming that 
the Universe follows the standard cosmology. Coherent oscillations begin in a radiation- 
dominated Universe, for which the Hubble parameter is 



i^(r) = /|^«.(T)-I1.1.66V54^-g-. (36) 

Here Mpi ~ 1.22 x 10^^ GeV is the Planck mass. We approximate the relativistic degrees 
of freedom 5'*(T) in the range of temperature we are interested in by 



UT) = { 



61.75, for T > A, 

10.75, for A > T > 4 MeV, (37) 
3.36, for T < 4 MeV. 



Eqs. ([5]), ([7]) and fl36l) provide the temperature Tf"^ at which the axion field begins to 
oscillate: 

ystd ^ 1 GeV g:"\Tr)a\ for Tf ^ < A, ^^^^ 
1 871 MeV (?;'^''(Tf for rf^ > A. 

For the entropy degrees of freedom g^s{T) we use the approximation 

(g,{T), forT>4MeV, 
9s{T) = I (39) 

13.91, forT<4MeV. 

Entropy conservation in the standard cosmology leads to the following relation between the 
scale factor a^^'^{T) and the temperature T: 

gl^g^ (T) T a'"^ (T) = constant. (40) 

The present axion number density nf'^{To) is computed from the number density at T^^'^, 
Eq. ([8]), assuming conservation of entropy and of the number of axions in a comoving volume. 



a^*^(To) J • ^^^^ 
The misalignment mechanism gives a contribution to the present axion energy density 
^std.mis _ man^a'^{TQ). In uuits of the critical density pcrit we have 



" ~ Pont g.s{Tr) ■ ^ ^ 
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Inserting the numerical values we obtain 



9.23 X 10-3 g:"\Tr) {Of /(^.)) fJil for fa > fa- 



7/6 



^std,mis^2 



for fa< fa, 



(43) 



Here h is the Hubble constant in units of 100 km s ^Mpc ^, and fa is the PQ scale at which 
the two expressions in Eq. (143|) are equal. For the values for g^(T) as in Eq. ([3] 



fa = 9.9 X 10^^ GeV 



String decays give a contribution to the present axion energy density 



(44) 



^std,str _ ^std ^: 



std ostd,mis 



0.164 n 



std, mis 



(45) 



where for a'^*'^ in Eq. ( |T6|) we have taken A^^ = 1, ^ 



std 



1, r 



^std 



0.8 



52| and C = 4.9, 



consistently with our previous work 6J]. In Section VI we discuss the modifications to 



47 



50 



53 



55 



77 



78| are used. 



Eq. (H5|) and to our results if the values f^^'^ = 70, ^^^'^ = 13 

The present axion energy density in the standard cosmology is then given by the sum of 
the misalignment mechanism and string decays contributions 



std 



std, mis 



+ n 



stdjStr 
a 



(46) 



If the PQ symmetry breaks after the end of inflation (Scenario I, fa < Hj/2it), there is 
only one PQ scale fa for which the totality of cold dark matter is made of axions. There 
correspondingly are also a single value for the axion mass rria from Eq. ([6]), and for the 
temperature Tf^'^, from Eq. (1381) . Using the observed value of flpmAh"^ in Eq. ([2]) and the 
expressions for fl^^'^ in Section II with {9f f{6i)) = 8.77, we have |64| 



ystd _ ^j_27 ± 0.25) X 10^" GeV 



m' 



std 



85 ± 3 /ieV, 



and 



Tf^ = 956MeV. 



(47) 
(48) 



(49) 



In Scenario II no axions from axionic strings are present. The parameter space is bounded 
by the non-detection of axion isocurvature fluctuations in the CMB spectrum, see Eq. ( |T2|) . 
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One parameter in Eq. (fT2|) can be eliminated by using the equahty between the axion energy 
density and the CDM energy density. In the standard cosmology, the isocurvature bound is 



2.89 X 10-5 ^/m)f!!i2, for fa < 9.9 X lO^e GeV, 



HlA2 < 



1.96 X 10-Va,(t 



(50) 



for 9.9 X 10^6 GeV < /„. 
It can be approximated by 

1.31 X 10-4 fa,i2, for fa < 7.6 X 10^0 GeV, 
Hi,i2 < { 2.89 X 10-5 /f/l2^ for 7.6 x 10^° GeV < fa < 9.9 x lO^^ GeV, 
1.96 X 10-4 for 9.9 x lO^^ GeV < fa, 



(51) 



where Hj 12 = i^//10^^GeV and fa^u = /a/10^^ GeV. The power-law dependence of the 
bound on fa^u changes twice: at fa ~ 10^^ GeV due to the effects of anharmonicities |64| . 
and at fa = 9.9 x 10^^ GeV, due to the change in the dependence of the axion mass with 
temperature. 



III. AXION CDM IN THE LOW TEMPERATURE REHEATING COSMOLOGY 
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In the low-temperature reheating (LTR) cosmology 



24, 



25 



26 



B,fla 



16 
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19 
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22 



27l |. the Universe after inflation is dominated by a massive decaying scalar 



field down to the reheating temperature Trh- The reheating temperature is defined as 



the temperature at w 



lich the decay width I\ of the scalar field is equal to the Hubble 



expansion rate H{T) 19 



21 



22 



23 



T, 



RH 



90 



1/4 



(52) 



At T > Trh, the Universe follows the LTR cosmology; at T = Trh, it transitions to 
the usual radiation-dominated era. At T < Trh, the Hubble expansion rate is given by 

T) depends on the scale factor a^^^{t) as in a 



Eq. fl36l) : at T > Trh, the Hubble rate H 
matter-dominated epoch 



17 



19 



21 



22|, 



H{T) = //(Trh) 



a^™(TRH) 



,LTR 



(T) 



3/2 



87r= 
"90 



RHy 



rp2 

fRH 
Mpi 



,LTR 



RH) 



,LTR 



(T) 



3/2 



for T > Trh- 
(53) 



Here H{Trb) is given by its expression in the radiation-dominated epoch. 
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Using the relation between the scale factor a^™(T) and the temperature T during the 
LTR epoch jl6|. 

^*(T)2/3 T^/^ a^™(T) = constant, for T > Trh, (54) 
and the usual relation during the radiation-dominated epoch, 

g,{Ty^^ T a'*^(T) = constant, for T < Trh, (55) 

we find the Hubble expansion rate 



90-9*iT)M]rr for T< Trh, 



H{T)= { V (5g) 



"49b.h{T).^, for T>Ti 



90 V- ; T^„Mpi 



RH, 



where 



9RHiT) = (57) 

The relation in Eq. (15^ shows that the evolution of the Universe during the LTR stage is 
non-adiabatic. 

In the standard cosmology, the axion field starts to oscillate at a temperature Tf'^ given 
by Eq. ([7]) with the standard expansion rate H{T) in the right-hand side. In the LTR 
cosmology, H{T) differs from the standard expression at T > Trh, and the axion field may 
start oscillating at a different temperature Tj'"™. 

More precisely, if the standard temperature Tf'^ is less than Trh, then the axion field 
starts to oscillate when the Universe is radiation-dominated. Moreover, since H{T) is the 
same in both cosmologies at T < Trh, the oscillations start at the temperature T^^^ = Tf*'^ 
given in Eq. ( l38l) . In this case the results of Section HTl apply. 

On the other hand, if Tf^'^ would be larger than Trh, then the LTR temperature Tj"'^^ 
will be smaller than T^^'^. In this case, the axion field starts to oscillate when the Universe 
is dominated by the decay of the massive scalar field 0. The temperature T^^"""^ follows from 
Eq. with H(T) given by the first line of Eq. ( l56i) . Since the dependence of H(T) on T 
steepens from T^ to T^ as T becomes greater than Trh, it follows that 

With TRH,MeV = Trh/McV, we find 

1/4 
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LTR 



Numerically, 

'351 MeVg-f{Tt^^) T^g^ev Cil\ for T^^ < A, 
160 MeV g-^TiTl^^'') T^gMeV fa}l\ for > A- 
To summarize, if Trh < A, 

'123 GeVg-'/\Tf^)f-ll\ 

351 MeV^-f (Tr^) T^Jg^eV /a^, 
[160 MeV,-r(Tr^)T^£Mev/MV 



(60) 



n 



LTR 



for Tf ^ < Trh, 
for Trh < TLT^ < A, 
for A < T^TR; 



(61) 



if Trh > A, 



LTR 



123 GeV^;^/^(Tf'i)/„"^/^ 



12 5 



871 MeVg:'''\Tf^)Ui' 



for Tfi < A, 

for A < Tf d < Trh, 



(62) 



160 MeV gnTiTl^"''') Ti^MeV Cii\ for Trh < Tf ^. 



Also the present axion energy density is modified from the standard case if Tf' > Trh- 
We examine the misalignment mechanism and the production in string decays separately. 
String decays give a contribution to the present axion energy density 



^LTR,str ^ ^LTR^LTR,mis ^ 0.032 fi^TR.mis 



(63) 



where we used the values Na = 1, f^™ = 0.27, = 16/27 and C = 4.9 in Eq. (fTel), 

consistently with the discussion in Section II B[ 

For the misalignment mechanism, the axion number density at the present time can be 
found from Eq. ([8]) with Ti = T^^^ using the conservation of axion number in a comoving 
volume, UaiT) oc a~^{T). This gives 



n^™(To) 



^-(^r) (^Slg^)' ' for < ^RH, 

forT->TRH. 



(64) 



Here na{Ti) is the function given in Eq. ([8]). One clearly has 

^(To) = nf (To) for Tf ^ < Trh. 



(65) 



For Tf*'^ > Trh, one obtains a different axion density. To understand the origin of the 
difference, it is convenient to introduce the ratio between the present density vi}^^(Tq) in 
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the LTR cosmology, and the present density n^*'^(To) in Eq. (HT!) that would ensue if the 
cosmology were standard at temperatures T > Trh- We write, for Tf'^ > Trh, 



where 



(66) 



is the standard-cosmology ratio of the comoving number of axions A^^"^^ at the temperature 
Ti^™ to the comoving number of axions A^'^*'^ at the temperature Tf**^, and 

is the ratio of the LTR-cosmology volume V^'^^ to the standard-cosmology volume V^^'^ at 
the temperature Tf''^^. 

The ratio N^"^^/N^^'^ accounts for the fact that coherent oscillations in the axion field 
start at a different temperature in the LTR cosmology compared to the standard cosmology. 
The ratio v^LTR^ystd accounts for the fact that at temperature T^"^^ the scale factors, and 
so the volumes, in the LTR and in the standard cosmologies differ due to entropy production 
from the decay of the scalar field in the LTR case. 

Using the relations between temperature and scale factor during the radiation and LTR 
epochs, Eqs. PUI) and fl5^ respectively, we find 

iv^™_ g.sirn fT^V ^gg^ 



and 

V^'^ g.siTKu) g^TF"") Ui"™; ■ ^ ^ 

From Ti^"^^ < Tf^ (see Eq. (EHD), we find that A^^™ > A^'^*'^. But for T^^™ > Trh, 

^LTR ^ ystd_ rjj_^g latter factor dominates, and n^™(To) is less than ?t,^*'^(To). 

The present axion energy density from the misalignment mechanism follows as, in units 

of the critical density, 

m n^^^(Tn) I for Tf ^ < Trh, 



Pcrit nstd.mis A^LTR yLTK 

"a ^std ystd ; lOr i >iRH- 



Here O^tfi'mis is given in Eq. (jH 
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The first line of Eq. (17T|) . valid for Tf'^ < Trh, is numerically equal to Eq. (H3|l . The 
second hne of Eq. (17T|) . valid for Tf'^ > Trh, is 

^LTR,mis^2 ^ 1 1-78 X 10-^ (0f /(^.)) g]i!\Tl^^'') fa% T|H,MeV, for < UTk^) . ^^^^ 
[7.46 X 10-8 {eU{ei)) fl,2TRHMeV, for /a > /a(TRH), 

where 

/a(TRH) = 5.69 X 10^^ GeVg^'fiT^^'') T^^^Mev (73) 

is the PQ scale at which the two lines in Eq. ( 172|) match. 

The present axion energy density in the LTR cosmology is given by the sum of the 
misalignment mechanism and the string decay contributions 

nstd.mis AtLTR yLTR , „,LTR\ f„„ T^std ^ rp 

"a N^td "tTStd- [i- + a J, lorJj^ > Jrh- 

Here, a^^'^ and a^™ are the values of the ratio pf^{TQ)/ p^^^{Tq) in Eq. ([TBI) in the standard 
and LTR cosmologies respectively. 




A. Results for LTR 

We now derive the regions of axion parameter space where the axion is 100% of the CDM 
in the LTR cosmology. We then compare them to the standard-cosmology regions. 

The CDM axion parameter space in the standard cosmology depends on the PQ energy 
scale fa (or alternatively the axion mass m^), the initial misalignment angle 6i and the 
Hubble parameter during inflation Hj. In the LTR cosmology an additional parameter is 
included, the reheating temperature Trh- 

If the PQ symmetry breaks after the end of inflation (Scenario I, fa < Hj/2tt), there is 
only one PQ scale fa for which the totality of cold dark matter is made of axions. There 
correspondingly is also a single value for the axion mass m^. In the LTR cosmology, using 
the observed value for flcDuh"^ in Eq. ([2]) and the expressions for fi^™ in this Section, we 
find 

/.^™ = (3.67 ± 0.11) X 10^^ GeV^i/l(T,L™) T'^/^.v, (75) 

and 

^ = 16.9 ± 0.5neV,-i/^(Tr^) T^gMeV (76) 
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FIG. 1: The Peccei-Quinn scale /a™ as a function of the reheating temperature Trh for the axion 
to be 100% of the CDM in Scenario I {fa < Hi/2tt). Also shown are the PQ scale /f"^ in the 
standard cosmology, and various constraints (shaded regions) . 

In Fig. □ we plot Z^^™ function of Trh- The jumps and kinks in the fY^^ line are 
due to the different values of 5'*(Trh) and g^ifTj"'^^) in Eq. fl57|) . There is also a (visually 
small) discontinuity between the /^"'"^ and /^**^ lines at Trh = Tt'^ = Tj""^^ due to different 
contributions from string decays. In fact, from Eqs. (143!) and (!74|) we have 



LTR 



1.32 g;'^'\Tn {ef f{e,)) {fT^f^Y^' (i + «^*'), 
1.32 ^."^/^^(Tr^) {9^ f{e,)) {f^jfyi^ ^ ^ (1 + 



for Tf" < Trh, 



a 



LTR\ 



for Tf ^ > Trh. 



(77) 

Equating the two lines in Eq. ([77]) at Trh = Tf^ = Tj^^^, where N^"^^ = N^^'^ and V"^"^^ = 
^std^ we obtain 

/ 1 _L ^std \ 6/7 

/."™(TRH = Tn=/r (^^^^j . (78) 

Inserting numerical values, /^"'"^(Trh = Tf**^) = 8.06 x 10^° GeV, which is slightly higher 
than /f ^. 
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In Fig. [T] we also shade out the following bounds: the bound from white dwarfs cooling 
times in Eq. f|T5l) : the indirect bound on fa from the non-detection of primordial gravitational 
waves arising from fa < Hi /2ti and Eq. (|T3l) (region labeled "Tensor Modes"); the bound 
on Trh from Big Bang Nucleosynthesis; and the bound from the ADMX experiment 80, 8l| 
excluding a KSVZ axion with a mass ma between 1.9 /xeV and 3.3 yueV. The dashed line 
marks the requirement that the axion starts to oscillate in the LTR cosmology, Trh < Tf"^, 
with Tf''^ given by Eq. (jUj). The ADMX bound can be rephrased as an exclusion bound 
for the reheating temperature Trh- Using the expression for the axion mass in Scenario 
I, Eq. ( 1761) . the ADMX result corresponds to an exclusion of the region 72 MeV < Trh < 
110 MeV, valid for KSVZ axions. 

Depending on Trh, fa^^ may differ from ff'^ in Eq. ( 1471) by orders of magnitude. The 
maximum value of fa^^ is achieved for Trh = 4 MeV and is, with (7*(Trh) = 10.75, /a™ = 
(8.58 ±0.25) X 10^^ GeV. This value is three orders of magnitude larger than f^^ in Eq. (H7!l . 
As discussed in Section |Vl these large values of fa correspond to axion masses that are beyond 
the reach of current DM axion search experiments. 

In Scenario II {fa > Hj/2tt), the parameter space is bounded by the non-detection of 
axion isocurvature fluctuations in the CMB spectrum, see Eq. f|T2l) . For Trh > Tf**^, the 
isocurvature bound has the same expression as in the standard cosmology, namely Eq. (l50l) . 
For Trh < Tf'^, we eliminate 9i in Eq. f|T2|) by equating ^cdm^^ with the expression for ^^^^ 
derived previously in this Section. The LTR isocurvature bound is then, for Trh < Tf***^, 



1.05 X 10-2 <7i'3(TL™) ^/mTnlueY fl%, for 9.2 x 10^^ T'^/^.v GeV < < /.(Trh) 

1-1/2 
RH.MeV' 



5.13 X 10-2 T-ji£,v, for AI^Rh) < fa- 



(79) 

Here A(Trh) = 5.69 x 10^^ GeV (?^1/'(TL™) T^j, ^^^^ is given by Eq. ([73]). The LTR isocur- 
vature bound can be approximated by 
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1.31 X 10-^ /„,i2, for fa < 9.2 x lO^^ T'^S^v GeV, 

1.05 X 10-2 gZiTY''') T-l^,^ fl% for 9.2 x 10^^ T-j^g.^ GeV < /„ < A(Trh), 

5.13 X 10-2 T^i/l,y, for fa > UTru), 

(80) 

As in the case of the standard cosmology, there are two changes in the power-law dependence 
of on fa,i2 in the LTR cosmology, the first one being at fa = 9.2 x 10^^ Tpj^^'^i^gv GeV 
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FIG. 2: In the LTR cosmology, axions are the 100% of the CDM in the white region on the left 
(limited by a different line for each Trh) and in the narrow bands marked by horizontal lines in 
the lower right triangle (one line for each Trh)- 

and the second one at fa = faiT^n)- Notice that at large /„, the isocurvature bound is 
independent of fa- 

When Trh = Tf'^, the LTR and the standard isocurvature bounds coincide. This happens 

for 

= 4.6 X 1021 T^l^^^ GeV, for Tf^ < A, 



= 5.3 X 1028 Tj 



-6 

RH,MeV 



(81) 



GeV, for Tf ^ > A 



Fig. [2] shows the regions of the parameter space {fa, Hj) where the axion is 100% of the 
cold dark matter in the LTR cosmology. The axion mass scale on the right is Eq. ([6]) with 

= 1. The region labeled "Tensor Modes" is excluded by the non-observation of tensor 
modes in the CMB fluctuations, Eq. (fT3l) . The region labeled "White Dwarfs Cooling Time" 
is excluded from astrophysical observations of white dwarfs cooling times for KSVZ axions, 
Eq. (fT5l) [73]. The line fa = Hj/27i divides the region where the PQ symmetry breaks after 
inflation (Scenario I, fa < Hi/27c) from the region where it breaks during inflation (Scenario 
II, fa > Hj/27r). 
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Trh Line in Fig. 1 

4 MeV (8.58 ± 0.25) x 10^^ Q^y 72 ±2 neV dotted line 

15MeV (2.64 ± 0.08) x lO^^ GeV 235 ± 7neV dot-dashed line 

150 MeV (1.23 ± 0.04) X 10^2 GeV 5.04 ± 0.15;ueV dashed line 

TABLE I: Values of f^^^ and m^"^^ for the axion to be 100% of CDM in the LTR cosmology. 

In the lower right region (Scenario I), the axion is the CDM particle if fa equals the value 
given by Eq. ( 1751) . Table I lists the values we plot. For comparison, we also plot the value 
for the standard cosmology, /f^ = 7.27 ± 0.25 x 10^° GeV, mf^ = 85 ± 3 /leV (thick line). 

In the upper left region (Scenario II), we plot the isocurvature bounds to the allowed 
parameter space for the standard cosmology (thick line) and for Trh = 4 MeV (dotted line), 
15 MeV (dot-dashed line) and 150 MeV (dashed line). For a given Trh, the isocurvature 
bound with the LTR cosmology lies below the standard line, because the entropy dilution 
term ~ (Trh/Ti)^ lowers the axion energy density. Thus, more parameter space is allowed 
for the axion to be 100% of the CDM in the LTR cosmology than in the standard cosmology. 

In the allowed region of parameter space for Scenario II, the axion can be 100% of the 
CDM provided the value of 9i is chosen appropriately. This value does not depend on H 



because in Scenario II ^ 6f. In the standard cosmology 6i is a function of fa only 64. |65| 



9 \ 6/7 

"cDMfe \ for f < f or 0- > 10-3 

n _ ) V 0.236 ; ' lOI ^ OI (7^ ^ lU , 

/«.12 - \ ) „ \ 2/3 . V«^J 

for fa > fa or < 10-3. 



0.0051 efj(e7) 

In the LTR cosmology, we find that 9i depends on both fa and Trh 

\ 2/3 

^ for fa < fa{Tnn) or 6^ > 17 T^'^' 



/«:12 - \ ^ „ X 1/2 . „ 



for fa > /a (Trh) or Oi < 17Tj^H,MeV 



7.46x10-8 e2/(0^)TRH,MoV, 

In Eq. dSSD we took ^.(T^™) = 61.75. 

We plot the relation between 9i and fa in Fig. [3] for the standard cosmology (thick line) 
and for Trh = 4 MeV (dotted line), 15 MeV (dot-dashed line) and 150 MeV (dashed line). 
As Trh decreases, one departs from the standard cosmology. The value of 9i at fixed fa, or 
of fa at fixed 9i, increases when Trh decreases. The largest departure occurs at the smallest 
value of Trh allowed by the BBN, Trh = 4 MeV. 
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FIG. 3: The initial misalignment angle 9i as a function of the Peccei-Quinn scale fa for the axion 
to be 100% of the CDM in Scenario II {fa > Hj/2tt): standard cosmology (black solid line), LTR 
cosmology with Trh = 4MeV (red dotted line), 15MeV (green dot-dashed line) or 150MeV (blue 
dashed line). 

In all these scenarios, the value of 6*, becomes uncomfortably small at large fa if one 
maintains that the initial angle 6i should not be tuned to a small value for the Peccei- 



Quinn mechanism to solve the strong CP problem [7l|. For the sake of illustration, if we 
decide that 6i = O.lvr is uncomfortably small, the PQ scale fa could comfortably be less 
than 5 x 10^^ GeV in the standard cosmology, but could comfortably be higher in the LTR 
cosmologies. For example, it could comfortably be as high as 5 x 10^^ GeV if Trh = 4 MeV. 
This relaxes the demand of 100% axion CDM on theoretical models that prefer fa of the 
order of the Grand Unification (GUT) scale ~ 10^^ GeV. 
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IV. AXION CDM IN THE KINATION COSMOLOGY 
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We now discuss axion cold dark matter in the kination cosmology [28|, l29|, l30|, l3l|, l32 



35 



36 



3d, 



38, 



391]. The Hubble parameter for this pre-BBN cosmology is (see Ref. 
H{T) = 



m) 



90 



■fl'kin(T) 



90 > Mpi ' 



forT > Tkin, 
forT < Tkin. 



(84) 



Here T^in is the temperature at which the Universe transitions from kination domination to 



radiation domination, and 



fl'kin(r) 



9l{T) 



(85) 



5'*(^km) 

Entropy is conserved during the kination cosmology, so the scale factor during kination 
a'*™(T) follows the same temperature dependence as a'^*'^(T), 



gl's{T)Ta''"\T) = constant. 



(86) 



For this reason, in this Section we do not make a distinction between the scale factor in the 
standard and in the kination cosmologies. We write a'^'°(T) = a^*'^(T) = a(T). 

If Tf^'^ < Tkin, coherent oscillations in the axion field start in the radiation-dominated 
Universe, at the temperature Tf'^ given in Eq. (!38l) . On the contrary, if Tf*'^ > Tkin, coherent 
oscillations in the axion field start when the Universe is in its kination stage. In this case, 
the temperature T^^^ at which axion oscillations begin is given by the following expression: 



2.48 GeV g^,l/\T^'n T^LmoV f^u, for Tf - < A 



T* 



kin 



-1/3 



331 MeV g^i;'\T^'') T^;!^Mey C12 . for > A. 



1/7 



(87) 



Here Tkin,MeV = Tkin/MeV. 

The axion energy density in the kination cosmology has contributions from string decays 
and from the misalignment mechanism. 

String decays give a contribution to the present axion energy density 



kin,str 



„ kin okin,mis 



0.42 f^'"n 



kin okin.mis 



(88) 



where to compute a in Eq. f[T6|) we used = 1, ^ = 2.06, C = 4.9 and r is given in 



Eq. (EHj) as 



^-'^--^0.81n(^) = — In (SB 



r 



2 
3 



1.6 

"3" ^H{T^^^ 



(89) 
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In the last expression we used the relation t oc 1/H{T), the fact that at the time of the PQ 
transition tpq the temperature of the Universe is T = /„, and the fact that at the time ti the 
corresponding temperature is Tj^™. Using the expression for the kination Hubble parameter 
in Eq. ( 18^ and neglecting the term ln(-i/ -^^jfM^) ~ 1, we obtain 



5kin(Tlkin) 

f'^-=1.6 1n(/,/T,'^-). (90) 

The temperature T^*^™ is greater than A for any value of Trh and any value of fa for which 
there are contributions from string decays (Scenario I, fa < Hi/2-n). Thus, using the 
expression for T]^™ in the second line of Eq. (IHTll . we obtain 
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= 57 + y In fa,l2 - lnTkin,MeV. (91) 

In the region of the parameters of interest for kination, f'^™ ~ 35. Thus, axions from strings 
dominate the total axion population, the energy density fi^i'^'^^tr being one order of magnitude 
larger than fi^™'™*^. We notice that this is opposite to what we obtained in the standard 
and LTR cosmologies, where the radiation of axions from axionic strings is a sub-dominant 
production mechanism for cold axions. 

The contribution from the misalignment mechanism results from the conservation of the 
axion number in a comoving volume, na{T) oc a~^(T). This gives 

-^(To)= )^{, (92) 

iMTn(^) , forrf^>rkin. 

Here na{Ti) is the function given in Eq. ([8]). One clearly has 

n^^^(To) = nf (To) for Tf <i < T^^. (93) 

For Tf''^ > Tkin, one obtains a different axion density. As for the LTR cosmology, we 
introduce the ratio between the present density n^™(To) in the kination cosmology, and the 
present density n^*'^(To) in Eq. pTl) that would ensue if the cosmology were standard at 
temperatures T > Tkm- We write, for Tf**^ > Tkm, 

<d(To) A^std ystd' ^ ) 

where N^'^'/N''^'^ and -l/kin^-j/std defined as follows. The ratio N^'^'/N''^'^ is the standard- 
cosmology ratio of the comoving number of axions N^^^ at the temperature T]^'° to the 
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comoving number of axions A^'^*'^ at the temperature Tf''^. Using Eq. (l86l) . we write it as 
^kin ^^(Tkm) /a(ri^'")V na{T^'^) g^T^'") fTf^^' 



The ratio v^km^^std ^j^g ratio of the kination-cosmology volume V"^™ to the standard- 
cosmology volume V^^*^ at the temperature T]^™, 

^ = (^a«td(2^kin) J = 1' (96) 
The last equality follows because no significant entropy is released during the kination stage 



33|, so a'^"(T) = a^*'i(T). 

The present axion energy density from the misalignment mechanism, in units of the 
critical density, is therefore 



. for Tf <i < Tkin, , ^ 



Ostd.mis 

^'a Tv^td ) ioi 1^ ^ kin- 



Inserting numerical values, the first line of Eq. (1971) is given by Eq. (H3|) . while the second 
line reads 

^kin,mis^2 ^ ^^5Q g-y\T^^:) {el fa,12 T^lueV (98) 

Due to the peculiar dependence of the Hubble rate with temperature in kination, H{T) ~ T^, 
there is no distinction in Eq. ([98]) between (]kin,mis jikin > ^ jikin < ^_ 

Finally, the present axion energy density in the kination cosmology is given by the sum 
of the misalignment mechanism and the string decay contributions 

ffifd.mis/^^ stdX for Tf'i < Tkin, 

j^^std,mis^(l^^kin)^ forTf^>Tki„. 

Here, a^*'^ and a"^'" are the values of the ratio p^*''(To)/p™'^(To) in Eq. (|T6l) in the standard 
and kination cosmologies respectively. 



A. Results for kination 



We now derive the regions of the axion parameter space where the axion is 100% of the 
CDM in the kination cosmology. We then compare them to the standard-cosmology regions. 
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The axion parameter space in kination cosmology depends on /„, Hj, 6i, and the addi- 
tional parameter Tkin- 

If the PQ symmetry breaks after the end of inflation (Scenario I, fa < Hj/2it), there is 
only one PQ scale for which the totality of cold dark matter is made of axions. There 
correspondingly is also a single value of the axion mass rria- In the kination cosmology, using 
the observed value of f2cDM^^ in Eq- and the expressions for fi^™ derived in this Section, 
we find 

= (7.9 ± 0.2) X 10^ GeV,y^(r,.) ..^lej^^^^r ' ^^''^ 

Ot + — m J^^^2 — in i kin.MeV 

and 

m^- = 739±22meV<7;i/2(Tkin)f'^'^r-iMev (101) 
In Eq. fllOOp we used the explicit expression for f'^™ derived in Section II B[ 

In Fig. m we plot as a function of Tkm- The function does not present jumps, 
because both (7*(T]^™) and 5'*(Tkin) do not change in the domain of /^™. The discontinuity 
between the and /^**^ lines at T^in = Tf'^ = T]^™ is due to different contributions from 
string decays. In fact, from Eqs. (H3ll and (l99ll we have 

= { -^^ ^ . ' - . " ^ . 

" ^kin "^kin 

"yystd" ystd 



1.32 g;'^'\Tn {ef fie,)) if:%y/^ (l + a^^'^), for < Th„, 

1.32 g;'/'\Tf-) {9f fm (/,k- ^ (1 + «'^^), for Tf ^ > Tki„. 



Equating the two lines in Eq. ffT02D at T^n = Tf^ = T]^'", where A^"^'"^ = A^"**^, we obtain 



l + a«*^^'/' 



/^"'^(Tkin = Tf <^) = ft' [yVo^^) ■ 

We find /„^'°(Tkin = Tf<^) = 2.04 x 10^ GeV. 

In Fig. m we also shade out the following bounds: the bound from white dwarfs cooling 
times in Eq. f|T5l) : the indirect bound on fa from the non-detection of primordial gravitational 
waves arising from fa < Hj /27c and Eq. (fT3l) (region labeled "Tensor Modes"); the bound on 
Tkin from Big Bang Nucleosynthesis; and the bound from the ADMX experiment excluding 
a KSVZ axion with a mass between 1.9 fieV and 3.3 fieV. The dashed line marks the 
requirement that the axion starts to oscillate in the kination cosmology, Tkin < Tf''^, with 
rpstd gjygj^ Eq. (m\. 

The PQ scale /^'° is orders of magnitude lower than the PQ scale ff^ in the standard 
cosmology. The low values of in comparison with ff^ is due to two different reasons. 
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FIG. 4: The Peccei-Quinn scale as a function of the kination temperature for the axion 
to be 100% of the CDM in Scenario I {fa < Hi/2tt). Also shown are the PQ scale /f"^ in the 
standard cosmology, and various constraints (shaded regions) . 

The first reason is that, since coherent oscillations of the axion field start later in the kination 
cosmology than in the standard cosmology, the initial comoving number of axions A^'^^'^ is 
higher than A^^*'^. The second reason is that the contribution from axionic strings to fij^™ in 
the kination cosmology is much higher than the same contribution to in the standard 
cosmology. Then, at a given PQ scale fa, the energy density nf'' > Jl^*^. A lower PQ scale 
is thus required in order to have the same CDM energy density ^cdm- 

The PQ scale /^'" can be so small as to violate the limit from the white dwarfs cooling 
time in Eq. f[T51) . This imposes the requirement T^in > 217 MeV if axions are 100% of the 
CDM. This requirement is more stringent than the BBN constraint T^in > 4 MeV. 

In Scenario II {fa > -ff//27r), the parameter space is bounded by the non-detection 
of axion isocurvature fluctuations in the CMB spectrum, Eq. flT2l) . For Tkin > Tf^, the 
isocurvature bound has the same expression, Eq. (!50|l . as in the standard cosmology. For 
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^kin < Tf^, we eliminate 9i in Eq. (fT2|) by equating r^cDM with the expression for 
derived in this Section. The resulting kination isocurvature bound for T^in < Tf''^ is 

Hj^u < 7.48 X 10-^ ^lW)fl% TlilMeW (104) 
This bound can be approximated by 

fl.31 X 10-Va,i2, for/„<3.26xl07GeVTki„,MeV, , , 

Hi,i2 = { (105) 

[ 7.48 X 10^7 fl% r,y^'MeV' for fa > 3.26 X 10^ GeV Tki„,MeV. 

Contrary to the cases of standard and LTR cosmologies, in the kination cosmology there 
is only one change in the power-law dependence of Hj^u on fa,i2, namely at /„ = 3.26 x 
10^ GeVTkin,MeV- This change is due to the effects of anharmonicities. 

When Tkin = Tf'^, the kination and the standard isocurvature bounds coincide. This 
happens for 

'fa = 4.6 X 10^1 T^l^^^ GeV, for < A, 

(1U6) 

= 5.3 X 1028 T^l^^^ GeV, for Tf ^ > A. 

Fig. [5] shows the regions of the parameter space {fa, Hj) where the axion is 100% of the 
cold dark matter in the kination cosmology. The axion mass scale on the right is Eq. ([6]) 
with = 1. The region labeled "Tensor Modes" is excluded by the non-observation of 
tensor modes in the CMB fluctuations, Eq. (fT3|l . The region labeled "White Dwarfs Cooling 
Time" is excluded from astrophysical observations of white dwarfs cooling times for KSVZ 



axions, Eq. (fT5|) . A similar bound from supernovae applies to other axion models [42(]. The 
line fa = Hi/2tt divides the region where the PQ symmetry breaks after inflation (Scenario 
I, fa < Hj/2tt) from the region where it breaks during inflation (Scenario II, fa > Hi/2ti). 

In the lower right region (Scenario I), the axion is the CDM particle if fa equals the value 
given by Eq. fllOOp . Table II lists the values we plot. Notice that the line at T^in = 4 MeV 
does not appear in this figure because it is excluded by the bound from white dwarfs cooling 
times. For comparison, we also plot the value for the standard cosmology, f^'^ = 7.27 ± 
0.25 X 10^° GeV, m'^'^ = 85 ± 3 fieV (thick line). 

In the upper left region (Scenario II), we plot the isocurvature bounds to the allowed 
parameter space for the standard cosmology (thick line) and for Tkin = 4 MeV (dotted 
line), 300 MeV (dot-dashed line) and 700 MeV (dashed line). For a given T^in, the kination 
isocurvature bound lies above the standard line. This is due to the fact that N^^"^ > N^^'^. 
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Line in Fig. [5] 



4MeV 
300 MeV 
700 MeV 



(3.8 ±0.1) X 10^ GeV 
(5.45 ± 0.2) X 10^ GeV 
(1.24 ±0.04) X 10^ GeV 

TABLE II: Values of and mjj'° for the axion to be 100% of CDM in the kination cosmology. 



I. 63 ± 0.05 eV 

II. 4±0.4meV 
5.0 ± 0.2meV 



dotted line 
dot-dashed line 
dashed line 



Thus, less parameter space is allowed for the axion to be 100% of the CDM in the kination 
cosmology than in the standard cosmology. 

In the allowed region of parameter space for Scenario II, the axion can be 100% of the 
CDM provided the value of 6i is chosen appropriately. This value does not depend on Hj, 
because in Scenario II ag ^ 6^. In the kination cosmology we find that 6i depends on both 
fa and Tkin as 

^CDM^^ ^'^^(^kin) ^kin,MeV 



fa, 



12 



1150 ef fie, 



(107) 
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FIG. 6: The initial misalignment angle 9i as a function of the Peccei-Quinn scale fa for the axion to 
be 100% of the CDM in Scenario II {fa > Hj/lir): standard cosmology (black solid line), kination 
cosmology with Tkin = 4MeV (red dotted line), 300MeV (green dot-dashed line) or 700MeV (blue 
dashed line). 

We plot this relation between fa and 6i in Fig. [6] for the standard cosmology (thick line) and 
for Tkin = 4 MeV (dotted line), 300 MeV (dot-dashed line) and 700 MeV (dashed line). As 
Tkin decreases, one departs from the standard cosmology. The value of 9i at fixed fa, or of 
fa at fixed 6i, decreases when Tkm decreases. This is opposite to the behavior in the LTR 
cosmology. The largest departure occurs at the smallest value of Tkm allowed by the BBN, 
Tkin = 4 MeV. 

The problem of 6i values uncomfortably small to solve the strong CP problem is more 
severe in the kination cosmology than in the standard cosmology. 
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V. DISCUSSION 



A. Comparison to previous work 



Axions in kination cosmology were studied only in Ref. 401] and only as hot dark matter 



(i.e. thermally produced in the hot primordial soup). To the extent of our knowledge, cold 



De: 



16 



ore. 



13, 



20 



22l |. but only to 



dark matter axions in kination cosmology were not examinee 

Axions in the LTR cosmology were studied before 
determine the cosmological bound on the PQ scale in what we call Scenario I, namely 
fa < Hj/2tt, in which the Peccei-Quinn symmetry breaks after the end of inflation. Our 
work can also be used to set an upper bound on the PQ scale by imposing 

Qa < ^CDM h\ (108) 

These bounds can be read off the figures and the equations in this paper, all of which 
represent the equation Vtah"^ < QcDAih"^- We remark that therefore our work extends 
previous papers in that we have examined also the region fa > Hj/27t, where the PQ 
symmetry breaks during inflation, have updated ficDM to the current observational value, 
have used an improved constraint on Trh from Big Bang Nucleosynthesis, and have included 
anharmonicities in the axion potential. 

Our numerical result for the highest allowed value of the PQ scale in the LTR cosmology 
/a™ — 8-58 X 10^^ GeV, obtained for Trh = 4 MeV, differs from previous authors. Stein- 
hardt and Turner showed that the entropy production due to the decay of a massive 
scalar field raises the maximum PQ scale to f^^^ ~ 10^*^ GeV, but they were corrected 



by Kawasaki et al. 20| for using the value of Ti in the standard cosmology instead of 



the LTR cosmology. Kawasaki et al. ^0|] used Ti in the LTR cosmology and obtained 
~ 10^^ GeV; however, they used VLah"^ = 1 and Trh = 1 MeV instead of the current 
value VtcDuh^ = 0.1131 ± 0.0034 and the current BBN bound Trh > 4 MeV. Giudice et 
al. find for the maximum PQ scale in the LTR cosmology the value Z^™ ~ 10^^ GeV, 
which is higher than ours for the same reasons as for Kawasaki et al. 
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B. Effects of changing the string decay parameters 

The computation of in Eq. (gT]), f^^^ in Eq. ([75]) and f^"" in Eq. (fTOOD strongly 
rehes on the model used to describe the axionic string evolution and the energy spectrum 
of emitted axions. 

In Section [TO] we discussed the dependence of f and ^ on the model for the axionic string 
oscillation and radiation spectrum. There we showed how these quantities in a modified 
cosmological scenario are related to their values in the standard cosmology. For the latter, 
we used the values f^*'^ = 0.8 and ^^^'^ = 1, obtained assuming that an axionic string radiates 



axions in a broad energy spectrum 



52l . I54J . I75l | and that the axionic string network is 



a global string network 52|, [5^, [Sj, |55|]. In the following we refer to these assumptions as 
Model A. With these values, the axionic string contribution to the total axion energy density 
in the standard and LTR cosmologies is subdominant compared to the contribution from 
the misalignment mechanism, while it is dominant in the kination cosmology. 

We now discuss the modification to the axion parameter space when we assume that 
axionic strings radiate axions in a narrow energy ^_syectrum 47|, l50|, l53|, l77|| , and that the 
axionic strings network is a local strings network [78]. In this case, f^*'^ = 70 and ^^^'^ = 13. 
We call this set of assumptions Model B. With these values, the contribution from strings 
to the axion energy density in the standard cosmology is dominant (A^^ = 1, ( = 4.9), 

j^std.str _ 200^f'^•'^^^ (109) 

This affects the value of 1)^*'^ = f^^t^.str _j_ ^std,mis^ which in Model B is about two hundred 
times higher than in Model A. As a consequence, with Model B, the value of the PQ scale 
jstd which n^^'^ = ficDM in Scenario I is 

/f d = 9 3 X iqS GeV. (110) 

This value is smaller than that in Eq. (147!) . obtained with Model A. It is of the order of the 
astrophysical constraint from white dwarfs cooling times in Eq. f|T5|) . 

We conclude that, depending on the model for the axionic string and its emission spec- 
trum, f^^'^ in the standard cosmology can range from the value (7.27 ± 0.25) x 10^° GeV in 
Eq. (1471) to the astrophysical bound from white dwarfs cooling times 4x10^ GeV (the errors 
come from the error on Qcdm only). Correspondingly, m^*'^ can range from 85 ± 3 ^ueV to 
15 ± 1 meV. 
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In non-standard cosmologies, when going to Model B, we must redo the calculation of r 
and ^ using the formulas in Section II B[ For the Model B value C,^^'^ = 13, Eq. flMl) gives 
c = (752 + l)/52 ~ 0.158 and Eq. ([32]) gives 

^ ~ 10 (^2 - 3/3 + 2^2.41/32 -3f3 + iy (for = 13). (Ill) 

The values of ^ in Model B are then C^^^ = 2.8 for the LTR cosmology {(3 = 2/3)and 
^kin _ ^YiQ kination cosmology {j3 = 1/3). 

For the parameter f, we turn to Eq. (!26|) . The results from Model A favor the fast- 
oscillating axionic strings model, which predicts f as given in Eq. fl28|) . Model B points 
toward a slow-oscillating axionic string for which f is given by Eq. (!27|) . Using Model B, 
for the illustrative case 6 = (10^^ GeV)^^ and Trh = Tkm = 4MeV, we obtain f ~ 70 in 
the standard cosmology, f'"™ ^ 20 in the LTR cosmology and f^^^ ~ 4300 in the kination 
cosmology. 

The LTR and kination axion energy densities from axionic strings in Model B are then, 
with A^d = 1 and C = 4.9, 

^LTR,str^2 _ H.S^LTR.mis ^]kin,str^2 _ 3 g ^ l]^^"''^^^ (112) 

The higher axionic string contributions in Model B with respect to Model A sensibly lower 
the values of the PQ scales fY^^ and for which the axion is 100% of the CDM. We have, 
taking Trh = Tkm = 4 MeV, 

/LTR = (143 ± 0.04) X 10^^ GeV and f^"" = (4.3 ± 0.1) x 10^ GeV. (113) 

C. Distinguishing non-standard cosmologies observationally 

Here we discuss how one might be able to distinguish different non-standard cosmologies 
before BBN using properties of the axion cold dark matter population. 

One may try to distinguish non-standard cosmologies by measuring both the axion CDM 
density ^cdm and the axion mass m^. However, one immediately runs into the following 
problem. 

Assume, for example, that the axion is found to be the main CDM component and the 
axion mass is measured at ~ 10^'^ eV. These facts can be ascribed to two different 
cosmological models. The first model involves the axion field evolving in the standard 
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cosmology, with the dominant contribution to the total axion energy density coming from 
axionic strings and only a tiny fraction from the misalignment mechanism, as Model B 
would predict. The second model involves a stage of kination before BBN lasting until 
^kin ~ 900 MeV, with the contribution from axionic strings and from the misalignment 
mechanism of the same order of magnitude, as in Model A. 

These uncertainties in the production of axion from strings decay prevent distinguishing 
non-standard cosmologies with this method alone. 

One may complement the measurements of ^cdm and with a measurement of the axion 
CDM velocity dispersion 5v. The latter allows non-standard cosmologies to be distinguished, 
at least in principle. The argument proceeds as follows. 

When axions start to oscillate at the temperature Ti, axions from vacuum realignment 
and axionic string decay have momentum dispersion of order the Hubble scale at Ti 



8p{T^)-H{T,). (114) 

The momentum dispersion scales with the scale factor as Sp{T) oc l/a(T). In the standard 
cosmology, the velocity dispersion at present is then 



In the kination cosmology, 

"'kin\ /„/'T^kin\\ 



It is clear that if one has measured m^, a measurement of 6v will give the value of T^in. 
Similarly, in the LTR cosmology. 



<zp) . f ^i!xy'\, X 10- ,„/s. (117) 

A difficulty in measuring 5v may arise from virialization of the axion population within 



galactic dark ha! 
space evolution 



OS, although it has been claimed that 5v would be preserved in the phase- 



VI. CONCLUSIONS 



In this paper we have examined the parameter regions in which the axion is 100% of 
the cold dark matter density in cosmologies that are non-standard before Big Bang nucle- 
osynthesis. We have recognized two ways in which these regions change in going from the 
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standard cosmology to the non-standard cases. If the Peccei-Quinn symmetry breaks after 
the end of inflation (Scenario I), the axion CDM regions shift to different values of the axion 
mass TTia (or of the corresponding PQ scale fa). If the PQ symmetry breaks during inflation 
(Scenario II), the axion CDM regions can shrink or expand according to the cosmological 
model. 

We have considered two different non-standard cosmologies that change the axion CDM 
regions in opposite directions. In the low temperature reheating (LTR) cosmology, the axion 
CDM regions shift to lower axion masses in Scenario I and expand in Scenario II. In the 
kination cosmology, the axion CDM regions shift to higher axion masses in Scenario I and 
shrink in Scenario II. 

Different axionic string models lead to different quantitative results, but the overall mod- 
ifications from the standard cosmology follow the same trend. 

We have also commented on the possibility to distinguish standard and non-standard 
cosmologies using observable properties of the axion CDM population. We have tentatively 
concluded that the axion velocity dispersion may be a good indicator of the cosmology before 
Big Bang nucleosynthesis. 
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